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Abstract 

The correlator of two field strengths is computed from an effective action for 
Yang Mills theories, which contains both gluons and an auxiliary antisymmetric 
tensor field for the field strength as local variables. This action allows to relate 
explicitly many different approaches to confinement, and is computed using 
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string theory proposed by Polyakov. 
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1 Introduction 



Many different methods have been used in order to characterize the infrared 
behaviour of Yang Mills theories: the Wilson loop [0J , the infrared behaviour of the 
gluon propagator 0, the monopole condensate from a dual Higgs theory || and 
corresponding analytic models ||, a confining string theory || and field strength 
correlators ||. The area law of the Wilson loop has been related to the infrared 
behaviour of the gluon propagator in , to monopole condensates in the dual Higgs 
theory in |||4]] and to field strength correlators in ||. More recently, field strength 
correlators have been studied in the dual Higgs theory in || and in the confining 
string theory in ||. 

In the present paper we continue the study of an effective action for Yang Mills 
theories, which contains both gluons and an auxiliary antisymmetric tensor field for 
the field strength as local variables It has the virtue of allowing for an explicit 
duality transformation of its abelian projection, and allows thus to relate explicitly 
a 1/g 4 behaviour of the gluon propagator to a monopole condensate in a dual Higgs 
theory. Moreover, the parameters of this effective action can be computed from the 
bare Yang Mills action by integrating the Wilsonian exact renormalization group 
equations. Note that antisymmetric tensor fields have already appeared frequently 
in effective descriptions of the infrared behaviour of Yang Mills theory, motivated, 
e.g., by confining string theories |J and, again, duality transformations of dual Higgs 
theories. 

The aim of the present paper is twofold: first, in chapter 2, we compute the 
correlator of two field strengths from the corresponding effective action. We find 
indeed an exponential decrease at large distances, in agreement with results on the 
lattice [11] , from the dual Higgs model || and from confining strings |§ . Notably we 
find that the parameters which characterize the slope of the exponential decrease are 
not the ones which parametrize a 1/g 4 behaviour of the gluon propagator. 

Second, we reconsider the Wilsonian renormalization group flow of the parameters 
of the Yang Mills action in the infrared regime in chapter 3. All previous approaches 
in this direction fl2| , |T3H have been plagued with the appearance of a Landau singu- 
larity of the running coupling constant, which did not allow the integration of the 
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renormalization group equations with respect to an infrared cutoff k 2 down to k 2 = 
(unless the running coupling is put to a finite value at k 2 = by hand). By taking a 
higher dimensional operator in the effective action (in the ghost sector) into account, 
we find that the Landau singularity is avoided. Moreover, we obtain an effective ac- 
tion in the deep infrared regime, which corresponds to the one of an effective abelian 
theory and resembles the action of a confining string theory ||. Conclusions and a 
discussion of the physical interpretation of this result will be given in chapter 4. 

2 Effective Action and the Field Strength Corre- 
lator 

In order to define the effective action with an auxiliary field for the field strength 
we start with the Yang Mills partition function 

e" G ( J *x) = 1 |p(A,c,c) e - s "'- s » +M «- c+ ^ . (2.1) 
Here Sym is the standard Yang Mills action, 

Sy M = \j 'd'xF; v F; v , (2.2) 



and S g is the gauge fixing and ghost part: 

1 /„ .„\2 



S g = J d A x [— (d,A;) 2 + d,c a (5 ac d, + g f abc 4) d 
Next we multiply the right-hand side of eq. (|2.1|) with 



(2.3) 



1 = 1. J VH e-i/^^-^) 2 . (2.4) 
In addition we add a source X" for the auxiliary field H® u and obtain 



e 

with 



-G(J,K, X ,X) = _L_ J Fj c> c) e -S(A,H)-S 3+ J-A + K-H + xc +X c ^ 



s(A, h) = J d'x Q (f; u ) 2 - 1 -f; v h; v + ~ (^) 2 ) . (2.6) 
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Note that one could perform the Gaussian integral over H in eq. (|2.5|) and obtain an 
equivalent formulation 

e -G(J,K, X ,x) = j_f C; Q e -S YM -S a+ J-A + K-F + K-K + xc +X c ( 2 .7) 

The effective action including the auxiliary field H is defined through the Legendre 
transform 

T(A,H } c,c)=G(J 1 K, X ,x) + J-A + K-H + x-c + x-c . (2.8) 

A general effective action contains an infinite number of terms with arbitrary powers 
in the fields and derivatives, which are just restricted by the need to satisfy the 
Slavnov- Taylor identities. 

Let us now consider an expansion of the A ^ and H^ u dependent part of V in powers 
of fields and derivatives: 



+ \ (D»H» V ) 2 + t {D , H, v f + -L (d. A,) 2 + . . . (2.9) 

Here H® is defined by H* = \s^v P( t H% a , and the covariant derivative D^, acting on 
fields ip a in the adjoint representation of the gauge group, by 

D^ a = d^ a + g f abc Al<p c ■ (2-10) 

In this approximation T(A,H) thus depends on 7 parameters Z, n, m, h, (3, g and 
a, which have to be computed within some non-perturbative scheme. In the next 
section we will discuss the Wilsonian exact renormalization group approach, but here 
we proceed by discussing several important properties of T(A, H). 

First, the "reduced" gauge coupling g appearing in the covariant derivative ( |2.10|) 
and in the non-abelian part of the field strength has no direct physical meaning. 
In order to define a physical gauge coupling one first has to eliminate if" from eq. 
Q2.9D by its equations of motion. One obtains 

T(A) = ^f- {F, v f + ... (2.11) 
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where the dots denote terms of higher order in the covariant derivatives which are 
induced by the terms ~ h, [3. Z e ff in ( 2.11] ) is given by 



n 2 



Z eff = Z-— . (2.12) 

m 2 

A physical gauge coupling (independent of field redefinitions of A^) is now given by 



Qphys — g/y Z e ff . (2.13) 

Next we study the two point functions (in momentum space) of the fields A a ^ and 
H®, as obtained from the effective action ( |2.9|) in the Landau gauge a — > 0. One 
finds 



5 2 r ^ / p tl Pv\ p * 2 

1 Oab V 9 ' P A\P 



6<Pi(-p)5<pj(p) J Aa Ah V P 



P ^) = 2 (7 P ? + f+ 7 R 4 > ( 2 - 14a ) 

p z (Zm z — n z ) + zpp 4 



6 2 T 



5<Pi(-p)5<Pj(p) . 



-iSab (Pp^pa ~ Pad up) • Pah{P 2 



■2 
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p z [Zm z — n z ) + Zpp 4 



5 2 r \ 1 



<U (5p« 5 aX - 5 P \ 5 aK ) Phh,i(p 2 



\5<Pi(-p)5<Pj(p) J H « H " 

+8ab {Sp K Pa P\ ~ 5pA Pa Pk + &a\ Pp Pk ~ S aK Pp P\) PhH,2 (p 2 ) 



hp 2 + m 2 



P (p 2 ) - Z(yk ~ ^ + n2 (2 14c) 

' (hp 2 + m 2 ) (p 2 (Zm 2 — n 2 ) + Zj3p 4 ) 



Let us first check the gluon propagator Pa in the case where T(A,H) of eq. Q2.9Q 
corresponds to the classical action S(A, H) of eq. ( |2.6| ). After rescaling H — > Aif in 
eq. ( |2.6| ) in order to give it the appropriate dimension of a bosonic field (with A equal 
to, e.g., an UV cutoff) this choice of T(A,H) corresponds to the following choice of 
its parameters Z, n, m, h, (3 and g: 
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Z = 2 , 

n = m = A , 

h = (3 = , 

3 = So • (2.15) 

Then one finds that P4Q9 2 ) is given by 1/p 2 (and Z e ff by 1) as it should be. 

On the other hand, Pa{p 2 ) enjoys a remarkable property if the three parameters 
Z, n and m are related such that 

9 

n 

Z eff = Z- — ^0 . (2.16) 
Then the gluon propagator Pa{p 2 ) becomes 

^ 2 ) = ^ , Al = f , (2.17) 
i.e. for p 2 << A 2 , the gluon propagator behaves like 1/p 4 . In [ID] we have emphasized 



that at the same time when eq. ( |2.16| ) is satisfied, the abelian projection of the action 



1|) allows for an explicit duality transformation, where the dual action corresponds 
to the one of an abelian Higgs model in the broken phase. The abelian projection 
corresponds simply to a vanishing of the reduced gauge coupling g. (In this case 
the effective action would be invariant under a new gauge symmetry of the form 
5A^ = A°, 5H^ U = d^A^i, were it not for the "gauge fixing term" (3/2(d IM H llI/ ) 2 in eq. 

(PD-) 

The fields appearing in the dual action are a dual abelian gauge field B^, a Gold- 
stone boson ip and a free massless scalar \ (whose origin can be traced back to the 
"gauge fixing term" ~ 0). The duality transformation is given by 



Z F, 



fJbV 



nH, 



Vh <9 M Ufa, , 



d v (p — nB u / V Zh 

d v X = \JP dp 

Assuming the relation fl2.16j ), the dual action T becomes 



(2.18) 
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f(B,*>) = j(f£) ; 



4 v '2 
The mass rh of the dual gauge field is given by 



(2.19) 



m 
~h 



(2.20) 



where we assumed again Q2.16Q to hold. 

Now we turn from previously published results to a new quantity of interest, the 
correlator of two field strengths F£ v . Due to the presence of a source K^ u for the 



field strength in the form ( 2.7 ) for the partition function this correlator can simply 
be expressed in terms of the functional G: 



W*) *J,(o)) 



5G 



5G 



S ab (5„ 5 V „ - V 5 vp ) 5 (x) + SK ^ {x) 5Kb ^ M 
5 2 G 



SKUx) 5KU0) 



(2.21) 



pa\"/ J K=J=0 



(Here we have omitted the Schwinger strings, which we assume not to significantly 
affect its x dependence.) 



After the Legendre transform (2.8) the correlator becomes in terms of the effective 
action V 



+ 



-5 ab (<W 8 vo - V 8 vp ) 5\x) + H« v (x) H%{Q) 
( 6 2 T V 1 



K 6<pi 8^J HUx)HbpA0) 
Assuming no vevs of the fields and A a , 



for 



5T 



ST 







(2.22) 



(2.23) 



the correlator is thus given by the Fourier transforms of the propagators given in 
eq. (|2.14cj ). Conventionally |6|,3,PII the correlator is decomposed into two Lorentz 
invariant functions D(x 2 ) and DAx 2 ): 
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(f;m KM) 



+ 



1 (_d_ 

2 [dx~, 



Sab (5^ 5 va - 8^ 5 vp ) D{x 2 
d 



{Xp $va ~ Xa S vp ) + — (x a 5p p - X p 5 p a ) 1 D x (x ) 



(2.24) 



After taking the rescaling of if by A into account, our expressions for the functions 
D(x 2 ) and Di(x 2 ) become after comparing ( [2.24j ) with (|2.14c|) : 



D{x 2 ) = A 2 / tSi eiPX P ^(P 2 ) - 



(2tt 

D ^- i ^!ik e " Ip ' , "^ ) (2 - 25) 



with Phh,i and Phh,2 as in eq. (|2.14c| ). Both functions D(x 2 ) and D\{x 2 ) have 
been measured on the lattice jpT] , and at least D(x 2 ) is well fitted by a decaying 
exponential 

D(x 2 ) ~ e~ M|:c| , M ~ 1 GeV , (2.26) 

for |x| < 1 fm. 

Subsequently we concentrate on the function D(x 2 ), whose non- vanishing is known 
to be responsible for the area law of the Wilson loop, if the cluster expansion converges 
||. Note that in our case, if T would correspond to the classical action S and hence 
h = 0, m 2 = A 2 , the propagator Phh,i is simply 1/A 2 and thus D(x 2 ) vanishes, since 
two 5-functions cancel each other in eq. ( |2.25|) . For a general action ( |2.9| ) D(x 2 ) is 
easily evaluated and one finds, for |x| ^ 0, 

where K\ is a Bessel function. Thus one obtains indeed an exponential decay (modulo 
powers) as in (|2.26|) , with 

i ji 

M = 7K ' (2 ' 28) 

The result ( |2.27|) agrees with the one obtained in [|] on the basis of a dual Abelian 
Higgs model, and in from a confining string theory. In all cases the inverse ex- 
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ponential decay length M, eq. (|2.28|) , coincides with the mass m of the dual gauge 
field, cf. the second equality in eq. ( |2.2(J| ). 

Since our effective action (|2~^ ) contains also the gluon field we can investigate, in 
how far the behaviour of D(x 2 ) is related to the gluon propagator. Surprisingly there 
is no direct relation: a 1/p 4 behaviour of the gluon propagator as in eq. ( p,17|) depends 
crucially on the relation (|2.16| ) to hold, i.e. on a relation between the parameters Z, 
n, and m. On the other hand, the result ( |2.27| ) does not at all depend on a relation 
like eq. (|2.16|) . Even if this relation holds, the parameter M in ( |2.28| ) is not directly 
related to the parameter A c characterizing the 1/p 4 behaviour of the gluon propagator, 
unless the parameters h and (3 in the effective action ( |2.9|) happen to be close to each 
other. (Actually, the function D\(x 2 ) does depend on the parameters Z, n and m. 



However, present lattice data [IT| do not yet allow to study this function in detail.) 
Thus we see that the different pictures of confinement - based on the gluon propagator 
or on the field strength correlator - are not directly related. 

Of course it can be argued that the result (|2.27|) is only a trivial consequence 
of our ansatz ( |2.9| ) for the effective action, whose parameters should eventually be 
computed from the Yang Mills Lagrangian. Note, however, that the action ( |2.9| ) 
implies the result ( |2.27| ) only if the parameters h and m 2 turn out to be non-zero, 
finite and positive. The result of computation of the parameters of the action within 
the context of the Wilsonian exact renormalization group approach will be presented 
in the next section. 



3 An Infrared Fixed Point of Exact Renormaliza- 
tion Group Equations for Yang Mills Theories 



In the recent years much progress has been made in applying exact renormaliza- 
tion group equations [I3J] to gauge theories P^ , P^|ff3H TT| . Also sources coupled to 
composite fields can be introduced in this formalism |1| which allows to apply it to 
the present case. 

The exact renormalization group approach requires the introduction of an "arti- 
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ficial" infrared cutoff k into the partition function Q2.1|) . Then one exploits the facts 
that the corresponding k dependent effective action T k becomes equal to the classical 
action 5* in the limit k — ► oo (up to additional terms determined by the modified 



Slavnov- Taylor identities [|T5| , p!6| ), and that an exact functional differential equation 
fixing the k dependence of can be derived. Integrating this Wilsonian exact renor- 
malization group equation from some large value k = A down to k = provides 
us with the physical effective action rfc = o in terms of the parameters of some "high 
energy" effective action I\. By construction no ultraviolet invergences appear in this 
approach, if both A and are assumed to be finite. 

To be concrete, in the present case with a source K coupled to the field strength 
F, one defines the functional G k ( J, K, x, X) including the infrared cutoff A; by a cor- 
responding modification of eq. (|2.7|): 



e -G k (J,K, x ,X) — _ J X>(A,C,C) e - s YM-S g -AS k +J-A+K-F+K-K+ x -c+x-c ^3 -Q 

where AS k implements the infrared cutoff for the gauge and ghost fields: 



AS k 



d 4 p 



O 2 ) K(p) + c a (-P) R k 9 (P 2 )c a (p) 



(3.2) 



The functions R k and Rg modify the gauge and ghost propagators such that modes 
with p 2 « k 2 are suppressed. Convenient choices are 



= if V + Q - l) P, P*) I Z e . p 2 /k 2 , 

e -p 2 /k 2 

R> 2 )=P \_ e - P y k * ■ (3-3) 

The effective action in the presence of the infrared cutoff k is again defined through 
the Legendre transform 

f k (A,H,c,c)=G k (J,K,x,x) + J-A + K-H + X -c + x-c . (3.4) 

Some expressions become more handsome when written in terms of T k which is given 
by with the infrared cutoff term ASk subtracted: 
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r k (A,H,c,c) = r k (A,H,c,c)-AS k . (3.5) 

From the path integral ( |3.1| ) and the Legendre transform ( |3.2| ) it is straightforward 
to derive the exact renormalization group equations jR-IM] 



^-Ai^ Rk{ H m-XJ ,, ' (3 ' 6) 

Here the fields (fi = (A a ^ H®, c a , c a ) denote all possible fields appearing as arguments 
of r fc or r fc , and the index i corresponds to the field type and Lorentz and gauge group 
indices. The matrix has non-vanishing matrix elements only in the subsectors 
(A^A") and (c a ,c a ). The inverse functional (8 2 Tk/ ' 8^pi8^p m )^i , however, has to be 
constructed on the complete space spanned by (A^, H^ v , c, c) including the auxiliary 
field H^. 

The right-hand side of eq. ( |3.6|) corresponds to a one loop diagram with an 
arbitrary number of vertices or external lines, and an insertion of d k Rz v or d k Rg into 
a gauge field or ghost propagator. The vertices and propagators have to be derived 
from the k dependent effective action IV For a given parametrization of r*. (or T^), 
non-linear differential equations for the k dependence of the corresponding parameters 
are obtained from eq. ( |3.6|) by comparing equal powers of fields and/or derivatives 
acting on the fields on both sides. 

Now we turn to a parametrization of the A^ and H^ u dependent part of T^. in 
the form of eq. (|2.9|) , where from now on the parameters Z, n, m, h and f3 depend 
on the infrared cutoff k. Differential equations describing the k dependence of these 
parameters are obtained from eq. (|3.6|) by considering terms quadratic in A^ or H^ u , 
and they have been derived in the Landau gauge in ref. ||10|| . 

On the right-hand side of these differential equations appears the "reduced" gauge 
coupling g, since all vertices (either from the non-abelian part of F* or from the 
covariant derivatives D^) are proportional to g. Thus the renormalization group 
equation describing the k dependence of g is also needed. 

In principle this renormalization group equation could be obtained from the three 
gluon vertex or the terms trilinear in A^ in eq. ( |3.6|) , but due to the large number of 
contributing diagrams it is much more convenient to introduce the ghost sector and 
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to consider the relations implied by the Slavnov- Taylor identities. 

The simplest non-trivial parametrization of the ghost sector of is given by 



Tf ost = Z g d,cD,c 



= Z g d^c a d^c a + Z g g f abc d^c a A b ^c c . (3.7) 

Note that the Slavnov- Taylor identities imply that the reduced coupling g in eq. 
( |3.7[ ) equals the coupling g implicit in eq. Q2.9p. (Here we neglect modifications of 
the Slavnov- Taylor identities for k 7^ [0,0 which vanish for k — > 0). Furthermore 



the ghost-gluon coupling receives no quantum corrections in the Landau gauge. In 
terms of the exact renormalization group equations this implies d k {Z g g) = or 

g- 1 d k g = -Z g l d k Z g . (3.8) 

Thus the k dependence of g can be obtained from the k dependence of the ghost 
wave function normalization Z gi which is much easier to compute. Within the present 
parametrization of T k the corresponding exact renormalization group equation has 



also been given in flip] , hence a closed set of differential equations describing the 
k dependence of the 6 parameters Z, n, m, h, (3 and g has been obtained. They 
have been integrated numerically, with boundary conditions such that at k = A 
T k (A, H, c, c) corresponds to the classical action (cf. eq. (|2.15|) ): 



Z(A) = 2 , n(A) = m(A) = A , 
h(A) = (3(A) = , g(A)=g , 

Z g (A) = 1 . (3.9) 

As a result of the renormalization group flow we found indeed that for small k Z e ff(k) 
(defined in terms of Z, n and m in eqs. ( p,12|) and (|2.16| )) vanishes. However, it was 



not possible to reach k = 0, since a Landau singularity appeared in g{k) at some 
small, but finite value of k = k c : within the present approximation of the right- 
hand side of the equation for d k g(k) is negative definite and proportional to g 2 (as in 
the case of one loop (3 function for Yang Mills theories), thus a Landau singularity 
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cannot be avoided. Clearly, this is an artifact of the neglect of the contributions 
of higher dimensional operators to the right-hand side of the exact renormalization 
group equations. 

Let us now go beyond the previous simple truncation of Tk and consider the effects 
of such higher dimensional operators. Since the running of g(k) is governed by the 
running of Z g (k), cf. eq. ( |3.8|) , we have to consider such contributions to the right- 
hand side of the equation for dkZ g (k). First, operators with the same powers of fields 
as in eq. ( |3.7| ), but involving higher derivatives, will not solve the problem: these 



have already been considered in the second of refs. H13] , and they do not modify 
the negative definiteness of dkg(k) (or the positive definiteness of dkZ g (k), implying 
Z g (k c ) = with k c finite). 

Thus we proceed by adding an operator of higher powers in the fields to the ghost 
sector of T k , which will contribute to the running of Z g (k). We replace Y 9 k host of eq. 
~ by 



T ghost = z^^D^c + Xd^cD^cF^F^) 
= Z g 8, c a d, c a + Z g g f abc d,c a Al c c 

+ \d,c a d,(c a F d p F d p )+\g f abc d,c a A\c c F d vp F d vp . (3.10) 

The action of the covariant derivative in the term ~ A in eq. ( |3.10| ) has been chosen 
such that the Slavnov- Taylor identities are still satisfied in a simple way: technically 
speaking, since the variation of Yf. host with respect to <9 M c a is still proportional to a 
total covariant derivative, the effective BRST variation of the gluon field A a ^ is still a 
total covariant derivative as before. 

In principle, also a term quartic in the ghost fields could have been added to 
Y^ host . However, the corresponding contribution to d^Z g involves only a diagram 
with an internal ghost propagator. This contribution is relatively small compared 
to the contribution obtained below in the limit where the gluon propagator becomes 
infrared singular, cf. our results below. 

After inserting Yf, host of eq. ( |3.10|) into the exact renormalization group equation 
( |3.6|) and expanding the right-hand side to the order d^c d^c one obtains the contri- 
butions to dkZ g which are shown diagrammatically in fig. 1. The "tadpole" diagram 
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involving a ghost-ghost-gluon-gluon vertex is proportional to the coupling A in eq. 



A priori a huge number of diagrams contribute to the running of d^X, alone 15 
to the order A . Since we are not (yet) interested in precise quantitative results, but 
rather in the essential features of the system of differential equations, we will only take 
the leading contributions into account: to the orders A 1 and A 2 (no higher orders in 
A exist) we consider only those diagrams, which are leading in the case of an infrared 
singular gluon propagator, i.e. for Z e ff — ► 0. 

Instead of including all the contributions to d^X to the order A , we will imitate 
these contributions by a small, but non-vanishing value of A at the starting point 
k = A. This procedure is justified, since one finds that for k — > A becomes extremely 
large, and all contributions ~ A to dkX become relatively negligeable. Also we have 
verified that the numerical results for k — > are practically independent of the starting 
point value of A provided it is small enough. This is just a manifestation of universality 
of the Wilsonian exact renormalization group flow, i.e. generally the results for k — ► 
depend only very weakly on the irrelevant couplings in Tk at the starting point k = A. 
Hence the diagrams which contribute to dkX are finally just those shown in fig. 2, 
which are of the orders A 1 and A 2 . 

In order to write down the resulting exact renormalization group equations for 
dkZ g and dkX we use the following notations: we need the gluon propagator function 
Pa(p 2 ) of eq. ( |2.14aj ), and in addition the ghost propagator Pc(p 2 ) given by 



P G (P 2 ) = ^ ■ (3.11) 

The presence of the infrared cutoff requires the replacement of Z by Z + R in eq. 
(g7L4|), and of Z g by Z g + R in eq. (gTg ) , with 



e -p 2 A 2 

R = l - e -P 2 /fc 2 ( 3 - 12 ) 

for the choices in eq. (|373|). go denotes the bare gauge coupling (appearing at the 
ghost gluon vertex), and, following the discussions below eq. (|3^ ) the coupling g 
equals go/Z g . For a SU(N) gauge group we then obtain from figs. 1 and 2: 
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dkZg 



q 2 dq 1 
16tt 2 



d k R{q 2 



3Ng 2 o 2 



PlPa 



P A P£) - 12N A g 4 Pi 



(3.13) 



q 2 dq 2 
16vr 2 



d k R(q 2 ) 



9Ng 



7,2 



\Z 2 q A P\ +4AV P\P, 



G 



(3.14) 



The presence of the term ~ A in Y 9 k host in eq. ( |3.1(J| ) affects also the exact renor- 
malization group equation for d k Z, but not the equations of the remaining parameters 
n, m, h and (3. For completeness we present these five equations below. Their right- 
hand sides receive contributions from the vertices from the non-abelian parts of F^ u , 
i.e. from the first two terms in T(A, H) of eq. fl2.9|) , and from the vertices due 
to the non-abelian parts of the covariant derivatives acting on H^ v or H®, which 
are thus proportional to h or (3. (These latter contributions are not very important 
numerically. ) 

Below Pa, Pah etc. denote the propagator functions of eqs. (2.14), with the 
replacements of Z by Z + R as below eq. ( |3.11| ). Because of the required development 
in powers of derivatives (or external momenta) the quantities p 2 dPA/dp 2 etc. appear 
frequently, and we define for convenience 



P'=p' 



dP 



P" 



V 



d 2 P 



(3.15) 



dp 2 ' (dp 2 ) 2 
for all propagator functions P A , Pah etc. With d k = k 2 d/dk 2 and for a SU(N) gauge 
group the five exact renormalization group equations are 



d, Z 



Ng 2 

167T 2 Jo 



dp 2 p A d k R k 



Z 2 [ H P 3 + wlp , A 



+n 2 (P A PA 



-uZPa^—PaPah + t^PahP'a + PahP'a + ~7^PaP'ah + PaP'ah , 

'ah \^P'ah + \p'ah) + Pah (f^4 + ^P'a + ~P'k 

+Pa (p~ 2 {jP'hh,i + i2 P hh,i) + ^Phh,2 + y^P'hh,2 + ^P'hh,2 
+Z(3P a PahP 2 (8P AH + 8P AH + 2P^) 
+n{3 (P 2 AH p 2 (2P AH - y P ah ~ P'ah) ~ PaPah {^P'hha + 
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-PaPahP 2 (3Phh,2 + fp' HH ,2 + P'hh^ ) 

—Prh,2 + 5P'hh,2 + P'hh,2^ ) 
+nhP A P AH (hP' HH ^ + \p'hh,) + h2p lHP 2 [Phh,i + + \Phh,x 



+PhP 2 AH p 2 [5P' HHA - "-P'hh+) + \Z] Pi 
Xfj2 dp 2 p 4 d k R k \ n P\(\z P A -\n P AH 



(3.16a) 



167T 2 JO 

+n(3P A (&p 2 P 2 AH + 3P A P HHA + ^p 2 P a Phh,2 



-nhP 2 A [~Phh,i + P'hh,i) ~ 8Z(3p 2 PiP AH 
-!3 2 P AH p 2 (^P a Phh,i + Zp 2 P 2 AH + Zp 2 PaPhh,2^J 



-h 2 p 2 P A P AH {^Phh,i + P'hh,i 
+(3hp 2 P A P AH (3P HHA + P' HH>1 ) 

^ J™ dp 2 p 4 d k R k [n 2 P\ - A(3np 2 P 2 A P AH - 3(0 + h)P 2 
+f3 2 P A (2p 4 P A P 2 AH +p 2 P A P HHA +p 4 P A P HHt2 ) 
+h 2 p 2 P 2 P HH , 1 ] , 

Ng 2 roo 

16^2 



(3.16b) 



J o dp 2 p 4 d k R k Z 2 P A P g -(P A + P g 



(3.16c) 
(3.16d) 



167T 2 Jo 



dp 2 p 4 d k R k [ n 2 P 2 p~ 2 (lp A + \P' A + 2 -P'{ 



+Pn( - P 2 (fp AH + f P' AH + \P'{ H ) - P A P AH (^P A + \p'i) ) 
+hnP A (2P A P AH + 2P A P'AH - 2P' A P AH ) 



+(3 2 (Pl + ^P'hh,i + o P 'hhA + ^ 2 ^ + \P'a + ^ 



+Pa PahP 2 (Ip'ah + \p1h) + P 2 aP 2 (fPHH,2 + 4P^, 2 + \P'I H ^ ) 

+/> 2 P 2 (?P M>1 + \p' hh , + ip^ 5l + ^ 2 P^, 2 ) 



+(3h( - P 2 a (3Phh,i + 2P^ il ) + P 2 AH P 2 (P A + 2P^) 
-P> 2 (5P^, 2 + 2P^ 



HH,2 



(3.16e) 



16tt 2 



f°^ 2 p 4 k R k \ n 2 P 2 A p~ 2 
Jo 1 



— Pa + -P'a + -K 
12 3 A 3 A 



16 



+Pn( - P\ {^P AH + \P' AH + 2 -P'i H ) - P A P AH (fa + 2 -P'i) ) 

+p\pi (Iphh A + \p' HHtl + \p'u) + py (fp A + \p' A + \pi) 

+py (j p *** + 2P '^ + \ p '^) + PaPa ^ 2 (t p ' AH + \ p 'ah) ) 

-fShP\ (3P HHA + 2P' HHA ) ] . (3.16f ) 

Note that all integrals are trivially ultraviolet finite since, with the present choice 
of R, dkR decreases exponentially for large p 2 , and infrared finiteness is ensured by 
the presence of the infrared cutoff terms in the propagators. The integrals in eqs. 
(3.16), and hence the integration of the exact renormalization group equations for all 
parameters Z, n, m, h, (3, Z g and A have to be performed numerically. 



It has already been verified in |T0|] that near the starting point k = A, with the 
parameters as in eq. (|3.9| ), the correct one loop (3 functions are obtained. Thus 
both Z e ff as defined in eq. ( p.!2[ ) and Z g decrease for decreasing k, i.e. the reduced 



coupling g = go/Zg increases. Also A increases for decreasing k, because the first term 



on the right-hand side of eq. (|3.13|) dominates. 

This trend continues for some orders of magnitude of decreasing k, during which 
Z e ff becomes extremely small. Then, at some critical scale k c , a violent transition 



occurs, which is entirely governed by the two renormalization group equations ( |3.13| ) 



and ( p. 14 ). (It is well known that coupled non-linear differential equations can lead 



to quasi- singular solutions.) First, the right-hand side of eq. ( [3.13j ) turns negative, 
since A continued to increase. As a consequence Z g increases for decreasing k, i.e. the 
reduced coupling g decreases. This decrease is very rapid; nearly instantly g drops to 
a value very close to 0. At some stage during this decrease, however, the second term 
on the right-hand side of eq. (|3.14j ) starts to dominate. From this point onwards also 



A decreases for decreasing k. 

Soon thereafter all evolution comes practically to an end, since all right-hand sides 
of the renormalization group equations become numerically tiny: the first terms are 
proportional to 1/Z g with Z g extremely large (remember that both g and Pq behave 
like 1/Z g ), and the second terms are suppressed by a power of Xk 4 for k — > 0. 

In figs. 3, 4 and 5 we show our results, obtained numerically with g(A) = g = 1.2 
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and A(A) = 0.01, for Z e ff(k), g{k) = g /Z g (k) and \(k) respectively. First, in fig. 3, 
we plot log(Z e ff) versus t = —log(k 2 /A 2 ). The evolution is from left to right: for 
k 2 = A 2 t is zero, and t increases for decreasing k. One recognizes the decrease of Z e ff 
up to k = k c ; for k below k c (or t above t c ) the evolution of Z e ff stops: the right-hand 
side of all renormalization group equations of Z, n, m, h and (3 are proportional to 
g 2 , and for t > t c g 2 is vanishingly small. 

This behaviour of g is shown in fig. 4, where we plot log(g) versus t. First, as 
described above, g increases until it becomes nearly singular. Then, at k — k c , it 
drops instantly to a tiny value and remains there subsequently. 

In fig. 5 we plot log(X) versus t. One observes its steady increase, its maximum 
at k = k c and its subsequent drop to a somewhat smaller value. 

We have checked the independence of our result on the initial values of g and A 
at the starting point k = A: we have varied g{A) = g between 1.0 and 1.5, and A(A) 
between 0.1 and 0.001 and obtained always qualitatively the same behaviour of g, A 
and the other parameters, just the scale k c varies correspondingly. 

We have not been able to describe the quasi- singular behaviour of g and A around 
the critical scale analytically. However, this behaviour becomes already apparent 
within a simplified toy model: it is possible to estimate the momentum integrals on 
the right-hand sides of eqs. (|3.13|) and ( |3.14|) , where the main contributions come 



from q 2 ~ k 2 ■ (—log(Z e ff)). Then these equations can approximately be written in 
terms of Z e jf, A and Z g (with g ~ 1/Z g , Pq ~ l/q 2 Z g ) as follows: 



a ry Cl C 2 A k 4 



Z eff Z g Z eff 

c 3 A 2 k 4 c 4 A 
d k A = — - — (3.17) 

Z eff 6 e ff Z g 

with Ci . . . C4 of 0(1). Already for constant Z e ff (note that Z e ff varies only weakly 
around k c ) the solutions for Z g and A of the simplified system (|3.17|) of coupled non- 



linear differential equations have the same properties as the ones shown in figs. 4 and 
5. This constitutes an independent check of our previous numerical results, and shows 
also that the qualitative features are independent of additional small contributions 
from neglected diagrams, i.e. the precise values of the constants c\ . . . C4. 

We can also insert the numerical values for Z e ff, Z g , A and k for k < k c into the 



18 



right-hand sides of the eqs. ( p. 17 ) and find that for all parameters P = {Z g , A} we 
have 



£^<10- • (3.18) 

The same relation holds for all the other parameters P = {Z,n,m, h, (3} for k < k c ; 
this explains, why there is no visible evolution in this regime. 

Strictly speaking, we have not yet obtained an analytic infrared fixed point of the 
exact renormalization group equations: the evolutions would vanish identically for 
k — > only for g = 0. However, with g ~ 10 -12 for k < k c we are so close to this fixed 
point that the tiny deviation plays no role for all practically purposes. (Still, the 
coupled system of exact renormalization group equations (|3.13|) , ( |3.14|) and (3.16) 
is so complicated that we did not manage to prove that analytically g{k) — > for 
k->0). 

It may be helpful to get some feeling for the critical scale k c where the sudden 



changes in the evolutions occur. In eq. ( 2.28 ) in section 2 we have introduced a scale 



M = m/y/h, which describes the exponential decay of the correlators of two field 
strengths F^ v and thus has some physical meaning. For the ratio k c /M we actually 
obtain, independently of the initial values of g and A, 

k c /M ~ KT 1 . (3.19) 
Let us now interpret our results in terms of the effective action T(A,H) in the 



form of eq. ( |2.9|) . We have to keep in mind that this truncation constitutes both 
a low energy approximation (in the sense that higher derivatives or higher powers 
of the momenta have been neglected) and a weak field approximation, since higher 
powers of and have been omitted. 

First, the fact that the reduced gauge coupling g practically vanishes turns the 
action into a free quadratic action of N abelian gauge fields A ^ and N antisymmetric 
tensor fields H^ u . This justifies a posteriori the abelian projection, which is required 
in order to make the duality transformation ( |2.18| ) feasible, and the omission of the 
Schwinger strings during the calculation of the correlator of two field strengths in 
eq. flUI). 

Second, we have introduced two physical dimensionful parameters in section 2: 
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A c = m/yf]3 characterizing the l/p 4 behaviour of the gluon propagator, and M = 
m/Vh characterizing the exponential decay of the correlator of two field strengths. 
Although the final results of the parameters m, h, and (3 for k < k c depend strongly 
on the initial value of the gauge coupling go, we always find h ~ /3 and hence A c ~ M, 
and a dependence of A c and M on the bare gauge coupling g and the starting point 
A (which plays implicitly the role of an ultraviolet cutoff) of the form 



k 2 ~ M 2 ~ A 2 e -^ 2 ' 11 ^ (l + 0(g 2 )) (3.20) 

as it should be. 

Next we consider the possible definition of a physical gauge coupling in eq. fl2.13| ), 



Qphys = 9/yZeff (describing, e.g., the low energy behaviour of a correlator of two 
static colored sources due to dressed one gluon exchange, cf. the second of refs. fHf ). 
Both g and Z e ff are tiny for k < k c , but to our surprise we find that also the ratio 
coresponding to g p h ys is very small (~ 10~ 4 ). This does not imply that Yang Mills 
theories are free theories, but rather that a momentum dependent gauge coupling 
a(q 2 ) tends to zero for q 2 —>■ (remember that we have neglected higher powers 
of derivatives in our action). Actually a similar result has been recently found on 



the lattice ||19|| , and it coincides again with the descriptions of confinement based on 
effective abelian models as the dual Higgs model in |4|], the "confining string theory" 
based on compact QED (with monopoles) in ||, and the gaussian approximation to 
field strength correlators in [H]. 



4 Discussion 

In this paper we have made progress in the description of Yang Mills theories 
at low energies in terms of an effective action, which contains both gluons and an 
auxiliary field for the field strength as local fields. It has been known before that, 
for Z e jf — > 0, this action describes a l/p 4 behaviour of the gluon propagator and, 
at the same time, its abelian projection allows for an explicit duality transformation 
relating it to a dual abelian Higgs model describing the condensation of monopoles. 

Here we have shown that it also allows for the computation of the correlator of 
two field strengths; the result (|2.27 ), however, has already been obtained before on 
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the basis of the dual Higgs model and the confining string theory ||. 

Furthermore we have reconsidered the computation of the corresponding effective 
action within the Wilsonian exact renormalization group approach. We have included 
a higher dimensional operator in the ghost sector, which is responsible for the running 
of the reduced gauge coupling g via the Slavnov- Taylor identities. (The considered 
operator is actually the only one to this order in F^ u consistent with the Slavnov- 
Taylor identities, if one requires that the BRST variation of the gluon field is still a 
total covariant derivative.) We find that now the Landau singularity in the running 
gauge coupling disappears; on the contrary, the running gauge coupling becomes tiny 
at small scales. Clearly, the sudden jump of the running coupling at the critical 
scale k c resembles a phase transition; it is not quite clear, however, whether the same 
phenomenon would also appear as a function of the temperature. 

Many properties of the resulting quasi-abelian effective action (always at low mo- 
menta, i.e. to lowest order in a derivative expansion) have already been discussed 
at the end of section 3. Let us add a final point, which is due to the smallness of 
Z eff = Z - n 2 /m 2 : for Z eff -> the first three terms ~ F 2 , F ■ H and H 2 in T(A, H) 
of eq. ( |2.9| ) become a perfect square {Z/A){F - mH/^Z) 2 . Then, due to the van- 
ishing of g, a field redefinition of the form H — > H' = H + \[ZF jm makes T(A, H') 
nearly independent of F or the gauge fields A: the term ~ h involving d^H^ u (for 
g = 0) is invariant under this redefinition of H due to the Bianchi identity, only in 
the term ~ f3 expressions of the form d^F^ v appear. Apart from this term ~ /3 the 
resulting action T(H') has then the simple form 

r(tf') = ^ (H» u f + ~ (0„ % v ) 2 (4.i) 

which coincides with the weak field limit of the action of the universal "confining 
string theory" in ||, where the antisymmetric tensor field couples to the surface 
bounded by the Wilson loop like the field strength. 

Clearly the relevance of the term ~ (3 has to be better understood; leaving it 
aside, a coherent picture emerges: in the action ( |4.1| ) is claimed to be universal 
in the sense that it does not depend on details (as the gauge group) of the confining 
gauge theory under consideration. Here we have obtained it from the integration of 
the Wilsonian exact renormalization group equations in the infrared regime, which 
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should also be described by universality classes. 

However, within the present approach it is straightforward to include higher 
derivative terms or momentum dependent couplings in the effective action (as, e.g., 
in the second of refs. ||13||); this allows to recover the perturbative behaviour of the 
Green functions at large momenta. The approach is thus quite unique in allowing 
for a coherent description of the effective action both in the perturbative and non- 
perturbative regimes. If our results, in particular the infrared limit in the form of 
the action (|4.1| ), can be confirmed by including higher dimensional operators in the 
effective action, a quantitative treatment of the non-perturbative regime within an 
expansion in powers of derivatives for a given Yang Mills theory seems to be within 
reach. 



Figure Captions 

Figure 1 : Diagrams contributing to the renormalization group flow dkZ g . Straight 
lines correspond to ghost propagators, twiggled lines to gluon propagators, and 
the crossed circle denotes an insertion of d^R with R as in eq. ( |3.12|) . 

Figure 2 : Diagrams contributing to the renormalization group flow dkX- 

Figure 3 : log(Z e ff) versus t = —log(k 2 /A 2 ). 

Figure 4 : £n(g) versus t = —log(k 2 /A 2 ). 

Figure 5 : in{\) versus t = —log(k 2 /A 2 ). 
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